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ABSTRACT 


Operator  valued  analogs  of  the  Chandrasekhar  H-function  that  occur  in  the 
study  of  neutron  transport  in  a slab  with  continuous  energy  dependence  and 
anisotropic  scattering  satisfy  a system  of  nonlinear  integral  equations.  An 
appropriate  Banach  space  setting  is  found  for  the  study  of  this  system.  We 
show  that  the  system  may  be  solved  by  iteration.  We  extend  the  domain  of 
analyticity  of  and  H^^  by  means  of  bifurcation  theory. 
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SIGNIFICANCE  AND  EXPLANATION 


i 

! 


The  Chandrasekhar  H-equation  is  the  following  non-linear  integral  equation: 

1 

H(g)  = 1 + ^H(u)  f -•  H ( v)  ( u)  d\)  , 

where  ip  is  given,  ^ is  a parameter,  and  H is  the  unknown.  This  equation 
is  of  importance  in  a variety  of  physical  situations  involving  transfer  of 
energy  by  radiative  scattering  processes,  e.g.  neutron  transport. 

The  H-equation  is  related  to  convolution  equations  of  the  following  type 

00 

(*)  f (x)  - / k(x  - y)f(y)dy  = g(x)  . 

0 

In  anisotropic  radiative  transfer,  the  following  generalization  of  (*) 
occurs 

00 

(**)  f(x,(D)  - / / k(x  - y,a),a)’ ) f (O)' )dw' dy  = g(x,a))  . 

0 Q 

The  equation  (**)  is  studied  in  this  paper  by  treating  it  as  an  equation 

of  the  form  (*) . We  must  then  interpret  k as  an  operator,  i.e. 

k(x  - y)f(y)  = / k(x  - y,co,ui' ) f (y,io' idu'  . 
fi 

An  explicit  formula  for  k is  given  by  equation  (1.8)  in  the  paper. 

We  consider  here  an  H-equation  in  which  H is  regarded  as  an  operator- 
valued function  on  L (f2)  , i.e. 

00 

H(ij)f  = f (ui)  - / h (u  ,u,(o ' ) f (o)' )da) ' . 

0 

This  generalization  appears  in  the  study  of  radiative  transfer  in  a slab  with 
anisotropic  scattering  and  continuous  energy  dependence.  We  show  that  we  may 
solve  our  equation  by  iteration  and  obtain  a bifurcation  result,  i.e.  we  show 
that  in  some  cases  there  are  two  solutions.  Numerical  and  asymptotic  results 
of  interest  in  transport  theory  may  follow  from  our  work. 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


OPERATOR-VALUED  CHANDRASEKHAR  H-FUNCTIONS 


C.  T.  Kelley 


I.  Introduction 

In  (10|,  Mulllkln  considered  Wiener-Hopf  factorizations  of  certain  operator-valued 
functions  that  arise  in  the  study  of  one-speed  neutron  transport  in  a slab  with  anisotropic 
scattering.  This  factorization  has  the  form 
(1.1)  (I  - = I . 

In  (1.1)  K is  the  Fourier  transform  of  a compact  operator-valued  function  K(x)  and 
I is  the  identity.  ~ ^ ~ ^ Fourier  transforms  of  compact  operator-valued 

functions  and  whose  supports  lie  in  (0,<»)  . ? is  a complex  parameter  and 

so  that 

(1-2)  sup  11k(X)  II  = 1 . 

X f ]R 

In  (1.2),  li  • II  denotes  spectral  radius.  Mulli)cin  derived  a system  of  coupled  non- 
sp 

linear  integral  equations  for  the  functions  H^(v,C)  and  Hjj(v,c),  where 
(1-3)  H .(v,c)  = H (i/v,c).  Rev  > 0 . 

r » *.  r I K,  ■“ 

The  purpose  of  this  paper  is  three-fold.  First  we  show  that  Mulliltin's  equations 
are  valid  with  only  minor  modifications  in  the  more  general  case  of  anisotopic  transport 
depending  continuously  on  energy.  We  develop  notation  that  allows  the  nonlinear  system 
for  H^  and  H^  to  be  written  in  a compact  way  in  a convenient  Banach  space  setting. 

We  then  show  that  the  equations  for  H^  and  may  be  solved  by  an  iterative  method. 

Mullilcin  and  the  author  [8]  have  recently  proved  a similar  result  in  the  case  of  isotropic 
scattering.  Earlier  results  of  this  type  may  be  found  in  [1],  (2),  and  [6].  The  proof 

given  here  differs  from  that  of  [8]  in  that  the  power  series  expansions  of  H^  and  H^ 

play  a direct  role. 

Finally,  we  show  that  H^  and  H^^  are  analytic  in  /l  - ^ for  C in  a cut 
neighborhood  of  1.  Our  proof  depends  on  a bifurcation  analysis  of  the  nonlinear  system 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024  and  the  National 
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for  and  Applications  of  results  of  this  nature  may  be  found  in  [7],  [9], 


and  [11]  . 

The  equation  for  steady  state  neutron  transport  in  a slab  with  anisotropic  scattering 
and  continuous  energy  dependence  may  be  written  as  follows,  112], 


(1.4) 


3iii 


V q—  (x,a))  + a;];(x,u))  = F(x,u) 
3x 


F(x,o))  = S(x,a))  + C / Kg  ) i|<  (x,(i)' )du)' 

n 


Here  ui  = [Eq'Ej^]  xS  =fi,  x > 0.  E denotes  energy  and  •p  a direction  vector  v is 


the  cosine  of  the  angle  made  by  the  vector  P and  the  positive  x-axis  a is  a continuous 


positive  function  of  E on  , o ^ 1.  K is  continuous  on  fl  x n and  satisfies 


(1.5)  The  integral  operator  on  C(S})  having  )cernel  Kg(u,(i))o  ^ (E ' ) 


is  positive  in  the  sense  of  [5]  and  has  spectral  radius  1. 
S(x,u))  is  continuous  on  [0,")  x n and  satisfies 


(1.6) 


/ sup  I S (x,u) |dx  < " . 

0 uen 


We  seek  solutions  ij/  that  satisfy  the  boundary  conditions 
(1.7)  iK0,a))  =0,  u ^ 0 


lim  i(i(x,u)  =0,  U £ 0 . 
x-«” 


We  let  C(Q)  and  L^(fl)  denote  the  spaces  of  continuous  and  essentially  bounded 
functions  on  0 with  the  sup-norm.  ]R  and  C are  the  real  and  complex  numbers. 

For  B a Banach  space  and  I an  interval  let  L(B)  be  the  algebra  of  bounded 


operators  on  B,  Com(fl)  the  algebra  of  compact  operators  on  B,  Bp(I,S)  (1  ^ p < “) 


the  space  of  Bochner  integrable  functions  from  I to  B. 


We  define  Ke  B,  ((O,*),  Com(L  (SI)))  by 
1 " 


(1.8) 


K(x)f(u))  = / K (w,a)')e"l’‘l°‘^’’^l^’ ' 0(xu')  -j^  f(w')dw' 

' ' 


In  (1.8),  0 denotes  the  Heavyside  function. 

We  may  then  convert  the  problem  given  by  (1-4)  and  (1.7)  to  a Wiener-Hopf  equation 


for  F € B,  ( (0,®)  , L (Q)  ) 
1 ® 


-2- 


(1.9)  F(x)  - C / K(x  - y)F(y)dY  = S(x),  x > 0 . 

0 

The  assumptions  on  and  the  theory  developed  in  [4],  (5)  and  [6]  imply  that  (1.9) 

has  a unique  solution  F for  |c|  < 1 and  that  the  function 
00 

(1  " CK(A))  = I - C / K(x)e^*^dx,  admits  the  Wiener-Hopf  factorization  (1.1).  Let 


S2^  = {(D  e n|vi  ^ O}.  We  state  some  properties  of  and  that  will  be  useful  in 

what  follows,  these  properties  follow  easily  from  [10]. 

Prop  (1.1).  Let  B(A,C)  be  any  of  the  operators  H (A,?),  H ^ (A,C),  then  the 

° « IT  ^ f IT 

following  hold 


(i) 

B(A,C)  - I 

is  an  analytic  Com(L^(fi))  valued  function  for 

|t|  <1,  ReA 

> 0 

(ii) 

B(A,C)  is 

a continuous  L (L_^(n) ) -valued 

function  for  | c | < 

1,  ReA  ^ 0 

(iii) 

B(A,t)  - I 

is  an  integral  operator  with 

)cernel  in  L (£5  x fi) 
00 

for  Id  < 

1 

ReA  ^ 0. 

Moreover  if 

is  given  by  (1.3),  the  )cernel 

h . (v,ui,u’ , (;)  of 
r,i 

Hr^/v.C)  - 

I 

is  nonnegative  on  f!  x fi 

, in  fact 

(iv) 

(v,<D,u',c)  > 0 on 

n X 

0 V ^ 0, 

0 < 

C < 1 

(V) 

h 

r 

(v,(o,u',?)  > 0 on 

fi  X 

V ^ 0, 

0 < 

c < 1 

(Vi) 

is  continuous  on 

[0,» 

) X n X n 

for 

Id  < 1 

(vii) 

h 

r 

is  continuous  on 

(O,” 

) X n X 

for 

Id  < 1 

(viii) 

h 

r 

and  are  increasing  functions 

of 

V for 

for  V ^ 0,  0 < ^ 

< 1. 

-3- 


r 


II. 


If  A is  an  integral  operator  with  kernel  in  ni  , we  denote  by  A’  that 

integral  operator  having  as  its  kernel  the  transpose  kernel  of  A.  We  define  I'  = I. 
We  let  N denote  the  algebra  of  integral  operators  on  Clfl)  having  kernels  that 


are  in  C(Q  x (2) . We  let  k.  (u,u')  be  the  kernel  of  A e N.  We  define 

A 


(2.1) 


'N 


(U 


st^ 

,u')^xn 


k (u,u ' ) 
A 


N is  a Banach  algebra  with  this  norm.  Note  that,  for  A,B  e N 

(2.2) 


k^(u,u’)  = / k^(w,a)Q)kg(Ujj,U))dajQ  . 


We  let  denote  the  algebra  of  integral  operators  on  b^((2)  having  kernels  which  are 


continuous  on  S2^  x f).  The  norm  on  N^  is  defined  by 


(2.3) 


)A|)  = sup  ) (k  (ii),u)’)  ) . 

+ u>en 

+ 


bi'ea 


As  in  [5] , [6]  and  [7]  we  consider  2x2  diagonal  matrices  of  the  form 


(2.4) 


A = 


° 

u 


A'  I 
r ' 


where  Aj^  and  A^  are  integral  operators  having  kernels  in  L^(R  x n)  . We  denote  the 


space  of  such  matrices  by  X and  define,  for  A e X 


(2.5) 

If  A,B  r X,  define  AB  and  A 


Alljj  = -"axdlk^  IKIlk^  II J . 

H r 


(2.6) 


AB  = 


A'B'  I 
r r ' 


' A 0 
r 


(2.7) 


A = 


Let  Xq  be  the  algebra  formed  by  adjoining  the  identity  to  X.  We  have 
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0 \ 0 I *0 


(2-9)  III  II  = 1 • 

^0  0 

When  no  confusion  can  result  we  drop  the  subscript  X_  and  refer  to  I as  I.  Let 
C([0,l),Xg)  = C be  the  space  of  continuous  X^  valued  functions.  The  norm  on  C is 


defined,  for  F e C,  by 


|f||  = sup  ||f(v)|L 

0<u<l  0 


Define  a map  L : C -*•  C{[0,1],X)  as  follows.  For  F(v)  e C write  F(v)  = a(v)I  + G(u) 
where  G(v)  e C([0,1),X),  then 

(2.11)  NlF)/^>  ° Kp(a..w-)e(-u’)  ; |^.|  a(|u'|) 

t / K^(a>,m^)e(-u„)  a(i^)vV|g  “'gJ  ^ “o'“i ‘'“o  ' 


k^I^).(v,(u,a,')  = a(u)K-(a),<o')e(U)  t |u| 


^ I a"(S)v  ^ y • 


Define  H(v,C)  = 


Hj^(v,5) 


H^(v,5) 


. We  have,  as  in  (10),  for  |?|  < 1,  Rev  > 0, 


(2.13)  H(v,C)  = I + tL(H  )(v,C)H(v,C)  . 

We  consider  this  equation  in  the  subspace  of  C given  by 

(2.14)  = {F(v)  = al  + G(v)  e C|Gj^  e n,  G^  e for  all  v,  a e T}  . 

We  may  then  consider  (2.13)  as  an  equation  in  the  space  C^.  We  define  T e X by 


The  assumptions  on  K imply  that  if  T is  considered  as  an  operator  on  C(!l)  x C(n) 


then  there  are  positive  functions  u , u.  in  C(fl)  so  that  the  vector  u = 

r i u. 


satisfies 


-5- 


u - Tu  = 0,  / u (u)u  (ai)d(u  = 1 . 


If  F ,F  e C((0,1J,X)  we  say  cil  + F > 61  + F if  a > 6 and  k > k almost  every- 
12  12  12 

where  in  (0,1  x SI  x n.  if  u = f 1 , v = f | e C(S!)  x (7(0)  we  say  u > v if  u > v , 

ri  I r£ ) “ 

^ everywhere  in  S!.  For  u e C(S1)  x C(Si)  and  a e K,  we  say  u ^ a if  u ^ |^j. 


Define  ff(C)  « by 


«(C)  = 


«^(0,C)  0 


W (0,0 

r 


Note  that 


ff  (0B(0  (I  - 5T)  =1,  Id  < 1 . 


The  following  lemma  is  crucial  to  what  follows.  The  proof  is  similar  to  that  in  (7) . 


Lemma  (2.1)  . The  limit 


lim  H(C)  = H(l) 


exists  in  C^.  Moreover,  H(i)  satisfies 

(2.19)  H(l)  = I + L(H*(1))H(1)  . 

Proof.  The  second  statement  clearly  follows  from  the  first. 

Define  a map  e 

(2.20)  LqF  = lim  (LF(v))  . 

The  definition  of  H and  (2.18)  imply 

//*‘^(C)  = I - C(Lq(H*))*  = ff(C)(I  - CT)  . 

For  0 < C <1  we  apply  both  sides  of  the  above  to  the  vector  u and  note  that  }i(0  ^ I 


to  get 


u - ;(Lq(H*)  )*u  1 (1  - C)u  . 


Set  M = - L,  we  have 

H*(v,C)(I  - C(Lq(H*))*)  = I - iH*<\/,0  (M(H*))*  . 


Hence,  as  H (v,^)  ^ 0, 


-6- 


u - H (v,?)(M(H  ))  u ' (1  - ?)H  (V,t)u 


H*(u,C)[(l  - C)U  + C(M(H*))*u]  < u . 


Let  V = (1  - C)u  + i;(M(H*))*u  = 


r 

‘'l  / 


. By  the  definition  of  M,  inf  > 0,  and 

ue  SJ 


inf  v^(ii))  - > 0,  therefore  h^  and  h are  bounded  on  S2  x f2  x [0,1)  uniformly  in  5 


As  0 < H(v,C)  £ H(v,Cq)  for  0 < C 1 Cq  < 1,  we 


have  that 


lim  h!  (v,u,u',C)  = h (v.ai.u')  exists  in  L ([0,11  ^ 0 x SI)  by  the  dominated 


convergence  theorem. 


To  complete  the  proof  one  may  now  proceed  exactly  as  in  [7] . 

The  lemma,  implies  that  H(?)  is  analytic  in  |c|  < 1 and  continuous 

on  Id  1 1- 

The  main  result  of  this  section  is  the  following 
Theorem  (2.1).  The  sequence  given  for  | C | £ 1 

(2.21)  ”0  “ ^ 

H ^ = I + CL{H*)H 

n+1  n n 

converges  to  H(C)  in  uniformly  for  |d  1 • 

Proof.  H(?)  is  the  unique  solution  of  (2.13)  which  is  analytic  in  l^l  < 1.  Let 

be  the  power  series  expansion  of  H(C)  about  C = 0.  We  have  that  i 0 fot 
n=0  ^ 

each  n,  P„  = I,  and,  for  n > 1,  P = I L(P.)P. . Let  S = I P . By  considering 

0 — n n . , K N " 

il+k*n  n*0 

the  coefficients  of  the  polynomials  , one  can  easily  show 

||h(c)  - h (c)  II  < I|h(UI)  - h (Id)  II  iI|h(i)  - H (1)  II  < ||h(i)  - s II  . 

n ““  n II  I* 

S converges  to  H(l)  by  lemma  (2.1),  and  the  fact  that  P„  1 0 for  each  n. 
n “ 


Hence,  the  proof  is  complete. 


III. 


In  this  section  we  show  that  H(v,C)  is  an  analytic  valued  function  of  *1  - C 


in  a cut  neighborhood  of  ; = 1.  We  do  this  by  showing  that  C = 1 is  a branch  point 
of  order  2 of  the  equation  (2.13).  The  proof  follows  the  general  theme  of  those  in 
[7],  (9],  and  [11]. 


In  this  section  set  H = H(l),  e = 1 - C,  and  G(e)  = (H  - H(c))H 
to  show  that  G satisfies,  for  t i 0 

e 


It  is  easy 


(3.1) 


G(e)  - HL((G(e)H)  ) = 


[(I  - G(e))H  - I)  - G(e)HL((G(£)H)  ) . 


1 - E 

The  Frechet  derivative  of  the  map  given  by  (3.1)  with  respect  to  G at  (G,e)  = (0,0) 


is  (I  - L),  where  L e given  by 


L(K)  = HL((KH)  ) , K e Cp  . 


(3.2) 

As  in  [7] , (9) , and  [11]  we  consider  an  alternate  operator  M defined  by 


(3.3) 


W(K)  = HM(  (KH)  ) , K e Cq  . 


In  (3.3)  M = Lq  - L.  We  require 


Lemma  (3.1) . There  is  P £ satisfying 


(i)  P = P 

(ii)  P is  independent  of  v 
(iii)  P > 0 


(iv)  P - A/(P)  = 0 

,2 


(V)  P = P . 

For  the  present  we  assume  the  lemma.  A computation  will  show  that  the  function  Q(v)  = vP 
satisfies  (I  - L)q  = 0.  Moreover  if  K is  such  that  (I  - W)K  = 0,  then  the  function 
S(u)  = vK  satisfies  (I  - L)S  = 0.  Bifurcation  of  (3.1)  at  E = 0 will  ta)te  place  if 
we  can  show  that  the  only  solutions  to  (I  - W)K  = 0 are  of  the  form  K = aP,  a £ I and 


that  the  range  of  {1  - M)  has  codimension  1 in  C^. 


For  K = 


K' 

r 


the  operator  M K has  the  form 


1 


\ 0 7 

2 2 

Moreover  and  (W^K^)  ' are  strictly  positive  elments  of  C([0,1],N)  and 

C([0,1J,N^)  respectively  if  Kj^  ^ 0,  0 and  ^0,  0.  As  is  compact  on 

C^,  we  have  that,  up  to  a constant  multiple,  P is  the  only  solution  to  (I  - M^)P  = 0 
and  hence  to  (I  - M)P  = 0 and  that  the  range  of  (I  - M)  has  codimension  1 in  C^. 
Therefore  Q(v)  = vP  is,  up  to  a constant  multiple,  the  unique  solution  to  (I  - L)Q  = 0. 
Moreover,  there  is  A « such  that  A((l  - L)K)  = 0 for  all  K e and  A(K)  > 0 if 

K ^ 0 and  K ^ 0.  One  may  then  proceed  as  in  [V]  , [9]  and  [111  to  obtain 
Theorem  (3.2).  H(i;)  is  analytic  in  /l  - 5 for  S sufficiently  near  1,  arg(l  - 5)  it. 


Moreover  we  have  the  expansion 


H(?)  = H - oe  ' QH  + 0(E)  , 


where  Q(v)  = vP,  and  P is  as  in  lemma  (3.1), 


(I  - L)Q  - 0 and  a = 


A(H  - I) 
A(vQ(v)) 


In  [7],  [9],  and  [11],  the  functional  A had  a simple  integral  representation.  In 
this  case,  however,  it  is  not  clear  that  A can  be  written  in  a simple  way. 

It  remains  to  prove  lemma  (3.1).  Let  u be  given  by  (2.16).  Consider  the  equality 

(3.6)  (1  - C)ff(0,i:)u  = (I  - C(Lg(H*(!;))*)u  . 

The  right  hand  side  of  (3.6)  has  a limit  as  5 approaches  1,  hence 

(3.7)  lim  (1  - C)ff(0,C)u  = w > 0 . 

C-^1 

kki 

* • ''r\ 

If  w = 0,  then  (I  - (L-(H  ))  )u  = 0 and  hence  there  is  v = , u 0 such  that 

0 V , 


♦ r 

V - Lg(H  )v  = 0 and  j v^  (u)  v^  (u) du  = 1 . 


Define,  for  x 


= T" 


e C(!))  X CiU) 


(3.9) 


Px  = 


f(x^,v^)Vr 


\(x^,Vr)v^ 


Then  P satisfies  (i)-(v).  If  u ^ 0 we  have,  by  (2.13), 

(1  - ?)ff(0,C)u  = (1  - C)u  + 5L|^(H*(C))  (1  - C)ff(0,5)u  . 

* 

Hence  w = L(H  )w.  Therefore  v satisfying  (3.8)  exists  and  P may  be  defined  by  (3.9). 

Theorem  (3.2)  implies,  in  fact,  that  w = 0,  and  hence  P is  unique  and  given  by 
(3.8)  and  (3.9)  . 
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